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Abstract 



The effective intercept parameter of the two-pion Bose-Einstein correlation 
function, A*, is found to be sensitive to the partial restoration of Ua(X) sym- 
metry in ultra-relativistic nuclear collisions. An increase in the yield of the 
vj meson, proposed earlier as a signal of partial Ua(1) restoration, is shown 
to create a "hole" in the low pt region of A*. A comparison with NA44 data 
from central S+Pb collisions at 200 AGeV is made and implications for future 
heavy ion experiments are discussed. 
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Introduction: Intensity interferometry is a useful method for studying the space-time 
geometry of high energy nucleus-nucleus collisions and elementary particle reactions (for re- 
cent reviews, see Ref. In particular, pion interferometry has proved useful in studying 
the space-time dependence of pion emission as was first shown experimentally by Gold- 
haber, Goldhaber, Lee and Pais ||. The method of intensity interferometry, known also as 
Hanbury-Brown-Twiss (HBT) correlations, was introduced by Hanbury-Brown and Twiss 
H for measuring the angular diameters of main sequence distant stars. The purpose of this 
Letter is to show that pion interferometry can be used to detect the axial {/a(1) restoration 
and the related increase of the rf production. 

As was shown in several papers ]5],|6| , at incident beam energies of 200 AGeV at the CERN 
SPS, the space-time structure of pion emission in high energy nucleus-nucleus collisions can 
be separated into two regions: the core and the halo. The pions which are emitted from the 
core or central region consist of two types. The first type is produced from a direct production 
mechanism such as the hadronization of wounded stringlike nucleons in the collision region. 
These pions rescatter as they flow outward with a rescattering time on the order of 1 fm/c. 
The second type is produced from the decays of short-lived hadronic resonances such as 
the p, N*, A and K*, whose decay time is also on the order of 1-2 fm/c. This core region 
is resolvable by Bose-Einstein correlation (BEC) measurements. The halo region, however, 
consists of the decay of long-lived hadronic resonances such as the u, rj, rj' and K s whose 
lifetime is greater than 20 fm/c. This halo region is not resolvable by BEC measurements. 
However, as will be summarized below, this region still affects the Bose-Einstein correlation 
function. 

In recent papers it was argued that the partial restoration of Ujl(1) symmetry of 

QCD and related decrease of the rf mass @-12| in regions of sufficiently hot and dense 



matter should manifest itself in the increased production of rf mesons. Estimates of Ref. 
show that the corresponding production cross section of the rf should be enhanced by a 
factor of 3 up to 50 relative to that for p + p collisions. 

The effective intercept parameter, A*, can be written in terms of the one-particle invari- 
ant momentum distributions j5|,|| of the core and halo pions and thus is sensitive to the 
abundance of the long-lived hadronic resonances such as the rf. To see this, consider the 
two-particle Bose-Einstein correlation function which is defined as; 

C(Ak,K)= A N ; iPl ^\ , (1) 
v ; iV 1 ( Pl )iV 1 (p 2 )' v ' 

where the inclusive n-particle invariant momentum distribution is given as 

1 do~ 

iVn(Pl,-,Pn) = —E 1 ...E n - — , (2) 

(T in (ipi...C{p n 

the relative and the mean four-momenta are given by 

Ak = p 1 -p 2 , K = — - — , (3) 

and p = (E p ,p). 

From the four assumptions made in the core-halo model fL3 |, the Bose-Einstein correla- 
tion function is found to be 
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C(Ak, K) ~ 1 + K(K = p; Q min )R c (Ak, K), 



(4) 



where the effective intercept parameter and the correlator of the core are defined, respec- 
tively, as 



and 



A* = A* (if = p; Q r 



R c (Ak,K) 



N c (p) 



N c (p) + N h (p) 
\S c (Ak,K)\ 2 



\S c (Ak = 0,K = p)\i 



(5) 



(6) 



Here, S c (Ak, K ) is the Fourier transform of the one-boson emission function, S c (x,p), and 
the subscripts c and h indicate the contributions from the core and the halo, respectively. 

In this form, \*(K = p,Q m in) is simply related to the momentum distributions of the 
core and halo pions. The Q m in dependence of A* which essentially indicates the separation 
of the core and the halo is actually defined by the experimental two-track resolution Q m i n . 

Axial Symmetry Restoration and the rj' production: In the chiral limit (m u = ma = m s = 
0), QCD possesses a U{3) chiral symmetry. When broken spontaneously, U(3) implies 
the existence of nine massless Goldstone bosons. In nature, however, there are only eight 
light pseudoscalar mesons, a discrepancy which is resolved by the Adler-Bell-Jackiw Ua{1) 
anomaly; the ninth would-be Goldstone boson gets a mass as a consequence of the nonzero 
density of topological charges in the QCD vacuum |T^Jl5| . In Refs. it is argued that the 
ninth ("prodigal" [^]) Goldstone boson, the rf, would be abundantly produced if sufficiently 
hot and dense hadronic matter is formed in nucleus-nucleus collisions. 

It was also observed, however, that the r)' decays are characterized by a small signal- 
to-background ratio in the direct two-photon decay mode. This may make the observation 
of rj' in this mode difficult, especially at small transverse momenta, where the increase is 
predicted to be the strongest. However, we now show that the momentum dependence of 
A* from pion correlations provides a good observable for partial Ua{X) restoration. 

If the v( mass is decreased, a large fraction of the rj's will not be able to leave the hot 
and dense region through thermal fluctuation since they need to compensate for the missing 
mass by large momentum |]?||8"|,|Tl~|| . These rj's will thus be trapped in the hot and dense region 
until it disappears, after which their mass becomes normal again and as a consequence of 
this mechanism, they will have small p t . The rj's then decay to pions via 



Tj — >• Tj + 



7T 



7T 



7T 



7T 



7T 



7T 



7T 



(7) 



Assuming a symmetric decay configuration (|pt|^+ — \pt\-K- — \pt\ij) an d letting = 958 
MeV, m v = 547 MeV and m n + = 140 MeV, the average p t of the pions from the rj' decay 
is found to be p t ~ 138 MeV. In the core-halo picture the r/', r\ decays contribute to the 
halo due to their large decay time (l/T^ > 20 fm/c). Thus, we expect a hole in the low 
p t region of the effective intercept parameter, A* = [N core (p)/N to tai(p)} 2 , centered around 
p t ~ 138 MeV. 

We note that the shape of A* will also be effected if the masses of the uj and rj mesons 
decrease. However, due to the large inelastic cross sections for u - meson scattering, the 
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lo are expected to rapidly reach chemical equilibrium when the hadronic fireball cools and 
their mass returns to its "normal" value. In this case, we do not expect a sizeable increase 
in the overall production of the u mesons. In addition, any enhanced production of the r\ 
mesons should only increase the depth of the hole primarily in the p t ~ 117 MeV region. In 
the case of equal production of the rj and //, there will be on the order of twice as many 7r + 
coming from the decay of the t/s than from the decay of the 77s. Thus, we concentrate on the 
dynamics of the rj' mesons giving an estimated lower bound on the depth of the produced 
hole. 

Description of the Simulation: In the following calculation of A*, we suppress the rapidity 
dependence by considering the central rapidity region, (—0.2 < y < 0.2). As a function of 



(8) 



m t = ypt+m 2 , A*(m t ) is given by 



K(m t ) 



N? tal( m t) 



where the numerator represents the invariant m t distribution of 7r + emitted from the core 
and where the denominator represents the invariant mt distribution of the total number of 
7r + emitted. The denominator may be explicitly written as 

= N :ore(mt) + N^ n+ (m t ) + AV^r+K) 

+AV^+ (mt) + N K o_ n+ (mt). (9) 

A detailed analysis |J has shown that the u does not to contribute to the core in the S+Pb 
reaction and in the NA44 acceptance. 

To calculate the n + contribution from the halo region, the bosons (u, rj', 77 and Kg) are 
given both a rapidity (—1.0 < y < 1.0) and an m t and then are decayed using Jetset 7.4 
16fl . The m t distribution 0,[nj of the bosons is given by 



N(m t ) = Cm?e- mt/T ^f, (10) 



where C is a normalization constant, where a = 1 — d/2 and where |T7|JT8| 

T eff = T fo + m(u t ) 2 . (11) 

In the above expression, d = 3 is the dimension of expansion, Tf Q = 140 MeV is the freeze-out 
temperature and (u t ) is the average transverse flow velocity. The m t distribution of the core 
pions is also obtained from Eqs. (TO) and flTTp. The contributions from the decay products 
of the different regions (halo and core) are then added together according to their respective 
fractions, allowing for the determination of \*{m t ). The respective fractions of pions are 



estimated from both the Fritiof |L9[ and the Relativistic Quantum Molecular Dynamics 
(RQMD) |2(| models as summarized in Ref. ||21|| . The calculation using Fritiof abundances 
is shown in Fig 2 (solid line). A similar m t dependence but with a slightly higher value of 
K(m t ) is obtained when using RQMD abundances. 

Simulating the presence of the hot and dense region involves including an additional 
relative fraction of 7/ with a medium modified p t spectrum. The p t spectrum of these 77' 
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is obtained by assuming energy conservation and zero longitudinal motion at the boundary 
between the two phases. This conservation of transverse mass at the boundary implies, 



*2 i *2 2 i 2 /in\ 

+ p w = m v , +Pt v >, (12) 



where the (*) denotes the rj' in the hot dense region. The p t distribution then becomes 
a twofold distribution. The first part of the distribution is from the rj' which have p* t < 



rn^i — m*?. These particles are given a p t = 0. The second part of the distribution comes 
from the rest of the r/'s which have big enough p t to leave the hot and dense region. These 
have the same, flow-motivated p t distribution as the other produced resonances and are 
given a pt according to the mt distribution 

N v r(m* t ) = Cm* t -°- 5 e- m * /r , (13) 

where C is a normalization constant and where T' = 200 MeV and m*, is the effective 
temperature and mass, respectively, of the hot and dense region. 

Assuming m*, = 500 MeV in the above scenario, the m t distribution of the 7r + from the 
decay of these rj' (rj' — > rj + 7r + + 7r _ ) is shown in Fig. 1. Also shown is the m t distribution 
of the 7r + from rj' assuming no hot and dense matter (Eqs. ( |I0"D and ( |TT| ) with Tf a = 140 
MeV, (ut) = 0.5 and m v ' = 958 MeV). Comparison of the two distributions shows the 
enhancement of the ir + in the low m t region which results from the presence of the hot and 
dense region. 

Using three different effective masses for the rf in the hot and dense region, calculations 
of A* (rat) including the hot and dense regions are compared to those assuming the standard 
abundances in Fig 2. The effective mass, ra* ; = 738 MeV, corresponds to an enhancement of 
the production cross section of the 77' by a factor of 3, while m*, = 403 MeV and m*, = 176 
MeV correspond to factors of 16 and 50, respectively. The two data points shown are taken 
from NA44 data on central S + Pb reactions at the CERN SPS with incident beam energy 



of 200 AGeV [22]]. The lowering of the 7]' mass and the partial chiral restoration result 
in a hole in the effective intercept parameter at low m t . This happens even for a modest 
enhancement of a factor of 3 in the rj' production. Similar results are obtained when using 
RQMD abundances. 

In addition, A*(m<) is calculated using Fritiof abundances with different average flow 
velocities in Fig 3. Here it is shown that A*(m t ) can also be a measure of the average 
collective flow. In our calculations, an average flow velocity of (u t ) = 0.50 results in an 
approximately flat, mj-independent shape for the effective intercept parameter A*(m f ), if 
the value of a = 1 — d/2 = —1/2 is kept fixed in Eq. (|TUp. Calculations using RQMD 
abundances result in a similar dependence on (u t ), but with slightly higher values of A* (rat). 

A limitation of our study is that we did not include the effects of possible partial coherence 
in the A* (rat) function. This is motivated by the success of completely chaotic Monte Carlo 
simulations in describing the measured two-particle correlation functions at the CERN SPS. 
However, a recent study P3| indicates that higher order BE symmetrization effects may 
also result in a decrease of A* (rat) at low p t . For the present system, this effect seems to 
be negligible, about a 1 % decrease, where the typical momentum scale of this effect is 



m t — m = T e ff and where the typical decrease is estimated JSJ from the measured radius 



and slope parameters. The flat shape of our A* (rat) distribution results from the inclusion of 
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the flow motivated temperature, T e ff, along with the effective, m t dependent volume factor 
PJl7|, V* oc m^ d ^ 2 } in Eq. (]T0|) . This flat shape reproduces the published NA44 data and 
differs from earlier theoretical calculations where A* is found only to increase with increasing 
m t . 

Summary: Our results reveal an important relationship between partial Ua(\) symmetry 
restoration and the shape (hole) of the A*(mt) parameter of the Bose-Einstein correlation 
function. We stress that this proposed signal is observed from the transverse mass de- 
pendence of the strength of the two-particle correlations, correlations which are presently 
being measured for fixed target Pb+Pb collisions at the CERN SPS. Measurements of two- 
particle correlations are also being planned for nuclear collisions at the Relativistic Heavy 
Ion Collider (RHIC) at BNL as well as at the CERN Large Hadron Collider (LHC). 

A qualitative analysis of NA44 S+Pb data suggests no visible sign of C/a(1) restoration 
at SPS energies. In addition, we deduce a mean transverse flow of (lit) ~ 0.50 in S+Pb 
reactions. Let us note that the suggested A*-hole signal of partial Ua(X) restoration cannot 
be faked in a conventional thermalized hadron gas scenario, as it is not possible to create 
significant fraction of the rj and if mesons with p t ~ in such a case. 
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FIGURES 



=L 70 




0.8 0.9 1 

m t (GeV) 

FIG. 1. Invariant m t distributions are shown for the ir + from the decay rf — > rj + vr + + 7r~. 
The solid line assumes the rj's come from a hot and dense region where the rj' have a twofold ^ 
distribution in which m*, = 500 MeV and T' = 200 MeV. The dashed line assumes the rj' come 
from nucleus-nucleus collisions at SPS energies where their mt spectrum is given by Eqs. (10) and 
(11), with T fo = 140 MeV, (u t ) = 0.5 and m v , = 958 MeV. 
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Effect of Partial U A (1) Restoration on X* 
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FIG. 2. Using the estimates of pion abundances given by Fritiof, the solid line represents 
A*(mt) assuming normal rf abundances while the other lines represent the inclusion of hot and 
dense regions, where T' = 200 MeV and the decreased mass of the rj' is m*, = 738 MeV (dashed 
line), m*, = 403 MeV (dotted line) and m*, = 176 MeV (dot-dashed line). These curves are 
calculated for (u t ) = 0.5. The data shown are from S+Pb reactions at 200 AGeV from the NA44 
collaboration. 
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<u t > Dependence of X* 
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FIG. 3. Using the estimates of pion abundances given by Fritiof, \*{mt) is calculated using 
different average flow velocities in the mt distribution. It is shown for (ut) = 0.00 by the solid line, 
for (u t ) = 0.25 by the dashed line, for (u t ) = 0.50 by the dotted line and for (u t ) = 0.75 by the 
dashed-dotted line. 



10 



